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In this lecture, we review the experimental situation of heavy Fermions with emphasis on the 
existence of a quantum phase transition (QPT) and related non-Fermi liquid (NFL) effects. We 
overview the Kondo lattice model (KLM) which is believed to describe the physics of those systems. 
After recalling the existing theories based on large-N expansion and various N=2 schemes, we present 
two alternative approaches: (i) a spin nuctuation-Kondo functional integral approach treating the 
spin-fluctuation and Kondo effects on an equal footing, and (ii) a supersymmetric theory enlarging 
the usual fermionic representation of the spin into a mixed fermionic-bosonic representation in order 
to describe the spin degrees of freedom as well as the Fermi-liquid type excitations. This kind of 
approaches may open up new prospects for the description of the critical phenomena associated to 
the quantum phase transition in Heavy-Fermion systems. 

Lecture given at the XXXVIII Cracow School of Theoretical Physics, Zakopane, June 1998: "New Quantum 
Phases, Elementary Excitations and Renormalization in High Energy and Condensed Matter Physics" (to appear 
in Acta Physica Polonica B 1999). 
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I. EXPERIMENTAL OVERVIEW 



In heavy- Fermion systems with 4f or 5f atoms (such as Ce or U) , the proximity of the electronic orbital to the 
Fermi level confers a Kondo effect at low temperature, i.e. an on-site compensation of localized magnetic moment 
by conduction electrons jl| . A direct consequence is the observation at low temperature of a very large effective 
electronic mass m* derived from the huge linear specific heat coefficient 7 = C /T and a correspondingly large 
Pauli susceptibility. Simultaneously, the realization of the de Haas- van Alphen (dH-vA) quantum oscillations [|| 
also concludes in favour of the existence of heavy quasiparticles. 

In addition to the Kondo effect, those systems are characterized by long-range RKKY (Ruderman-Kittcl- 
Kasuya-Yoshida) interactions between neighboring local moments mediated by conduction electrons. The com- 
petition between the Kondo effect and the RKKY interactions leads to the possibility of either a non-magnetic 
or a long-range magnetically-ordered ground state ||. A zero-temperature quantum phase transition occurs 
governed by the value of the exchange coupling J between the spin of the conduction electron and the local 
moment. One of the most striking properties of the heavy- Fermion compounds discovered these last years is 
the experimental possiblity to explore this quantum phase transition by varying the composition change 
(as in CeCuQ- x Au x or Ce x La\- x Ru-iSi-i), or by applying a pressure or a magnetic field. Thus a magnetic 
instability is observed at x c = 0.1 in CeCu§- x Au x Q, and x c = 0.08 in Ce x Lai^ x Ru2Si2 0- For x = x c where 
Tn = 0, the observed behavior at low temperature is at odds with that usually expected for a simple Fermi liquid 
(FL). In CeCue- x Au x , the specific heat C depends on T as C/T ~ —Ln(T/T ), the magnetic susceptibility as 
X ~ 1 — ctyT, and the T-dependent part of the resistivity as Ap ~ T (instead of C/T ~ ^ ~ Const and Ap ~ T 2 
in the Fermi liquid state). Pressure or large magnetic fields are found to restore the FL behavior. 
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The origin of this non- Fermi liquid (NFL) regime is a largely discussed problem. The three main interpretations 
which have been proposed rely on (i) a single impurity multichannel Kondo effect in which the internal degree 
of freedom is provided by the 4f or 5f quadrupolar moment |1, (ii) a distribution of the Kondo coupling due to 
the disorder leading to a distribution of the Kondo temperature P(Tk) @ and (iii) the proximity of a quantum 
phase transition pO|- p^[ as is emphasized in this course. 

Another important insight is provided by the Inelastic Neutron Scattering (INS) experiments carried out in 
systems close to the magnetic instability. The measurements performed in pure compounds CeCug or CeRu2S%2 
]l6||r^ ] have shown the presence of two distinct contributions to the dynamic magnetic structure factor: a q- 
independent quasielastic component, and a strongly q-dependent inelastic contribution peaked at the value cj max 
of the frequency The same experiments carried out in systems with varying concentrations as Ce x La\- x RuiSii 
[ p"8[ , show a shift of c<j max to zero when getting near the magnetic instability. Any theory aimed to describe 
the quantum critical phenomena in heavy-Fermion compounds should account for the so-quoted behavior of the 
dynamical spin susceptibility. 



II. THEORETICAL OVERVIEW 



The model which is believed to describe the heavy-Fermion systems is the Periodic Anderson Model (PAM) 
defined for the case of spin 1/2 as 

h = E e * c t c ^ + E "J2 flf* + v J2( c l^ + c - c -) + u T, n \i n fi (!) 

ha ia ia i 

where — f^fia- It describes the conduction electrons c^a with dispersion Sk which hybridize with the 
localized / electrons of energy Eq. The grand canonical ensemble is used and both energies Sk and E$ are 
measured from the chemical potential p. The hybridization matrix element is approximated by a constant V. U 
represents the on-site Coulomb repulsion between / electrons. 

In the Kondo limit, a canonical transformation allows to change the Periodic Anderson Model into the Kondo 
lattice model (KLM) defined as 

F = E£fc4 CT C feCT + jJ2 S i' S i ( 2 ) 

ka i 

where J is the on-site Kondo coupling between the spin of the conduction electrons Si — ^2 c la T aa'Cia> an d 

aa' 

the localized spin represented by Sj = j\ a T aa' ]ia' hi the Abrikosov pseudo-fermionic representation of the 

aa' 

spin imposing the constraint riji — = 1. That canonical transformation analogous to the Schrieffer-Wolf 

transformation for the single-impurity case is valid in the regime \ek\ << (—Eg) and \ek\ << \Eq + U\. One gets: 

J = V 2 — -g- + E +u . In the U — > oo limit, one has J — —V 2 /Eq while in the symmetric Anderson model 

defined by Eq = —U/2, the result is J = AV 2 /U. The Kondo lattice Model has been first introduced by Doniach 
in 1977 H and we refer to the paper of Tsunetsugu et al for an extensive review on the KLM essentially at 
D=l. 

Let us first recall the main physical ideas behind the Kondo Lattice model: (?) the competition between the 
Kondo effect and the RKKY interactions leading to the so-called "Doniach phase diagram" and (ii) the nature 
of the screening of the local moments in the lattice. 

Concerning the point (i), the competition between the Kondo effect on each site which tends to suppress 
the magnetic moment with decreasing temperature and the RKKY interactions which, on the contrary, tend 
to magnetically order the local moments, leads to the well-known Doniach phase diagram ||. Let us call Tjf 
the Kondo temperature for the single- impurity and (or Tq) the Neel (or Curie) temperatures in the absence 
of Kondo effect: Tg- = Z?exp (— l/p n J) and ~ {pqJ) 2 where D and p are the bandwidth and the density 
of states at the Fermi level of the conduction electron band. Thus, at small poJ, Tjy is larger than and a 
long-range magnetic order is established with eventually a reduction of the magnetic moment due to the Kondo 
effect. Oppositely, at large poJ, T% is larger than T^, the Kondo effect wins and the system does not order 
magnetically. Therefore, the real Neel temperature Tn first increases with increasing poJ, passes through a 
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maximum and finally goes to zero at a critical value of the coupling poJc giving rise to a zero-temperature 
quantum critical point. 

As far as the nature of the screening is concerned (point (ii)), an important idea advanced by Nozieres [ pl| 
is the possibility of an exhaustion of the conduction electrons in the screening of the local moments. In the 
single-impurity Kondo case, at low temperature, the local moment is screened by the conduction electrons and a 
spin-singlet state is formed. In this so-called Kondo effect, the number of conduction electrons in the screening 
cloud formed around the impurity is equal to 1. In the Kondo lattice case, the number of conduction electrons 
which are available are to be taken within a thermal window of width ksT around the Fermi level. It should be 
compared to the number of sites N$ to be screened. Depending on the value of the parameters, some situations 
may occur where the available conduction electrons are "exhausted" before achieving complete screening leaving 
residual unscreened spin degrees of freedom on the impurities. That idea of "uncomplete" Kondo effect is at the 
root of the supersymmetric theory that we propose later on. 

At high temperature, pertubation techniques in J may be applied leading to the famous Kondo minimum in 
the resistivity as a function of temperature. Those pertubation techniques fail below the Kondo temperature 
and there is a need for other techniques to solve the problem at low-temperature. The other approaches based 
on the Bethe Ansatz and the Renormalizaion Group which revealed very powerful in the single-impurity case 
cannot be generalized to the lattice case. In that context, there has been an intensive search for new approaches 
among which the large-N and more generally the functional integral approaches described here. 



A. The Large-N expansion 



An important breakthrough in the understanding of the periodic Anderson Model occured about 15 years ago 
when the idea of slave-bosons was introduced: |2^,|23| for the single-impurity case, 24 27) for the lattice. In the 



limit of large on-site Coulomb repulsion (U —> oo) where the double-occupancy is energetically forbidden, one 
can introduce a slave-boson representation in which the two allowed states i.e. the empty or the singly-occupied 
states are represented by ej|0 > and fj a \0 >. The exclusion of double-occupancy is expressed as 

Pi = ete, + -1 = (3) 

The physical electron creation operator which creates transitions between empty and singly-occupied sites is 
represented by / icr ej, while the number of / electrons on site i of spin a is equal to f irj eie\fi a = f irj fi a provided 
that the local constraint is satisfied. Hence the slave-boson representation of the U — ► oo PAM hamiltonian is 
given by 

h = E e * c L<*«r + E oJ2 flf* + v £(4«4/* + c - c -) ( 4 ) 

k , CJ i , (7 i , <7 

provided that the local constraint Eq.(^) is satisfied that is enforced with the aid of a time-independent 
Lagrange multiplier A^. The operators e\ and obey bosonic and fermionic statistics respectively. 

It is then convenient to generalize the original PAM model from SU(2) to SU(N) by allowing for the spin 
index a in Eq.(^) to run from — S to S. That corresponds to the situation of impurities of spin S coupled to 
conduction electrons of degeneracy N with N = 25+1 . The corresponding SU(N) model is interesting because 
it can be solved exactly in the limit N — > oo. For this limit to make sense, the hybridization matrix element V 
should be scaled as 1/yN therefore V = V/yN. 

We do not give here all the details of the calculations which can be found in the litterature. Let us say 
that the saddle-point approximation which consists to take and Aj as site-independent (as well as time- 
independent for a) is exact in the limit N — > oo.lt leads to the formation of two quasiparticle bands of energies 

= \ek + Sf ± \J (ejfe — £/) 2 + W 2 e^ where Sf = E + X . The values of eo, £/ and fi are fixed by the saddle- 
point equations. The two bands are split by a hybridization gap and the density of states at the Fermi level 
is strongly renormalized p(Ep) ~ 1/Tr- where Tr- = £>exp(£'o/po^ 2 )- The specific heat coefficient 7 and the 
magnetic susceptibility \ are a ls° strongly enhanced with a Wilson ratio x/l equal to 1. Oppositely, the charge 
susceptibility is found to be unenhanced. The ground state corresponds to a collective Kondo screening in which 
only a fraction of conduction electons equal to Tk/D screens each of the impurities. 
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Next step is to include the gaussian fluctuations around the saddle-point. The corresponding corrections in 1/N 
generate effective interactions among the quasiparticles which can be analyzed in terms of Landau paramaters. 
In the case of the multichannel single-Kondo impurity problem, non-crossing approximation methods |2j| have 
been extensively used to derive the finite temperature behavior with very accurate comparison with the Bethe- 
Ansatz and Conformal Field theory results. The generalization to the case of the lattice has still to be done. All 
the results can be reproduced by starting instead from the N — > oo Kondo lattice hamiltonian and performing 
a Hubbard Stratonovich transformation on the coupling term making the field $ appear. There is then a one- 
to-one equivalence between <£>o = JJ2ka < c \a-fka > and Veo. In all cases, no magnetic instability is found at 
the order 1 /N since the RKKY interactions occur at the order 1 /N 2 [^9| . The latter point constitutes a serious 
drawback of the large- N expansion which makes it inappropriate to describe the Quantum Critical Point observed 
experimentally. The fact that the slave-boson does not carry spin implies that spin and charge fluctuations 
cannot be treated on an equal footing contrary to what happens with other slave-boson representations as the 
one introduced by Kotliar and Ruckenstein that we present now. 



B. The different N=2 approaches 



In the case of a degeneracy N = 2 (S = 1/2), Kotliar and Ruckenstein (KR) |}0| introduced 4 slave-bosons 
Bi, pi a and di in order to keep track of the 4 possible local configurations so that the empty \0 the singly- 
occupied \a >i and the doubly-occupied | Ti>i states are represented by: \0 >i— e\\vac >, \a >,= p icr fl a \vac > 
and | Ti>i = d\f^f^\vac > where e;, pi a and di are bosons and /j CT fermions. Then the PAM can be written as 



(5) 



fe.CT 



with 



= (4put +Pl-a d i) I \j 1 ~ 4 e * ~ p\-aP*-<?\J 1 - d i d i - PiaPi* 



provided that the 3 following constraints are fulfilled 

Pi = et ei +^pJ^ CT -l = 



(6) 



Qia = fiJia - (PiaPia + d\di) =0 

The choice of the denominator in guarantees to recover the free electron gaz limit at U — > 0. This 
representation first introduced in the case of the Hubbard model was shown at the saddle-point level to give 
back the variational Gutzwiller approximation (GA) as developed by Rice and Ueda |n| . It then allows to 
include the gaussian fluctuations around the GA solution Q. In the case of the PAM @||,|9|^|], the KR 
representation already leads to interesting results at the saddle-point level as soon as staggered symmetry-broken 
state appropriate for bipartite lattice with nesting is allowed. Notably, at d=I |Q, the approach essentially gives 
the same results as those obtained by the variational wave function approach of Gulacsi, Strack and Vollhardt J35| . 
At infinite dimension |3fj] , a phase transition to the antiferromagnetic insulator is found below a critical value Vc 
of the hybridization matrix element consistent with the Doniach predictions and in quantitative agreement with 
the d = oo QMC J3lj and exact diagonalization |3f| results. The general phase diagram for the three-dimensional 
case has been determined by Doradzinski and Spalek f39f . Those studies enlighten on the nature of the moment 
compensation which takes place in the antiferromagnetic state. The study of the V-dependence of the staggered 
magnetizations m/ and m c shows an almost total moment compensation of and m c near Vc suggesting an 
itinerant magnetism in which f and c electrons are part of the same quasiparticles. Oppositely, in the V — > 
limit, rrif saturates while m c goes to zero indicating a local moment magnetism. The figure 4 of the paper [ |39[ 
illustrates the latter point by showing the V-dependence of m/ and m c in the antiferromagnetic insulating state. 
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Finally, other path integral approaches to the KLM have been proposed based on different Hubbard 
Stratonovich decouplings of the exchange term J4l| , [42|] . For the extended KLM in which the Heisenberg in- 
teractions among neighboring sites are included, we will mention the work of Coleman and Andrei [[l2| which 
consists to keep the Resonant- Valence-Bond (RVB) parameter Xij — JJ2aflafj<? 011 neighboring sites at the 
same time as the Kondo parameter $ = JJ2a c \afi<? quoted before. This approach leads to the stabilization of 
a spin-liquid state at low temperature with possible anisotropic superconducting instability. This approach has 
been used in a recent paper by Iglesias, Lacroix and Coqblin p3[ | where they propose a revisited version of the 
Doniach phase diagram in which the Kondo temperature is drastically reduced resulting of the formation of the 
resonant valence bonds. 

In the rest of the paper, we will develop two alternative approaches to the Kondo lattice model: (i) the first one 
consists to keep the f and c magnetizations at the same level as the Kondo parameter Jl5|] . We will show how it is 
possible to account for the spin-fluctuation and the Kondo effects on an equal footing thus combining both large 
N and spin-fluctuation theories. To our point of view, this approach constitutes an ideal framework to study 
the quantum critical phenomena around the magnetic transition; (ii) the second approach consists in enlarging 
the usual Abrikosov pseudo-fermionic representation of the spin into a mixed fermionic-bosonic representation in 
order to describe the spin degrees of freedom as well as the Fermi-liquid type excitations 44 1. The analogy of the 
approach with the supersymmetry theory of disordered systems leads to give it the nickname of " supersymmctric 
approach" . 



III. THE SPIN FLUCTUATION-KONDO FUNCTIONAL INTEGRAL APPROACH 



In the grand canonical ensemble, the hamiltonian of the Kondo lattice model (KLM) constituted by a periodic 
array of Kondo impurities with an average number of conduction electrons per site n c is written as 



ka 



(7) 



ka 



in which r are the Pauli matrices (t x ,t v ,t z ) and r the unit matrix; J is the antiferromagnetic Kondo 
interaction (J > 0). 

We use the Abrikosov pseudo-fermion representation of the spin S^: Si = fl a T aa i fi a > . The projection into 

<7<7 ! 

the physical subspace is implemented by a local constraint 

ia 

A Lagrange multiplier Aj is introduced to enforce the local constraint Qi. Since [Qi,H] = 0, \ is time- 
independent. 

In this representation, the partition function of the KLM can be expressed as a functional integral over the 
coherent states of the fermion fields 



VdaVficrdXi exp 



C{r)dT 



(9) 



where the Lagrangian C(t) is given by 

£(r) = £ (t) + H (t) + Hj(t) 



/ j 1(7 

ia 



H (r) 



£ fc c L c ^ - v N s 

ka 



V S ka ) i 
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Hj(r) = J^S/i-Sc 



With S Ci = c \a T <?<?' c i<y' and S /i = S i 
aa' 

We perform a Hubbard-Stratonovich transformation on the Kondo interaction term Hj(t). Since more than 
one field is implied in the transformation, an uncertainty is left on the way of decoupling. We propose to remove 
it in the following way. First, we note that Hj(t) may also be written as 



(10) 



where S/ Ci = ^2 ]\ a Taa'C%o' and nj Ci = ^ j\ a T® a ,c ia i (respectively S c / ; and n c f i their hermitian conjugate). 



The Kondo interaction term is then given by any linear combination of J ^2 S/i ' S ci (with a weighting factor 

i 

x) and of the term appearing in the right-hand side of Equation ( |To| ) (with a weighting factor (1-x)). x is chosen 
so as to recover the usual results obtained within the slave-boson theories pi|-p7t . One can check that this is the 
case for x = 1/3. The Kondo interaction term is then given by 



Hj(t) = J s J2( S f* ■ S c* + S /c, • S c /J - Jc^n fCt n cfz 



(11) 



with Js = J/4 and Jc = J/3. 

Performing a generalized Hubbard-Stratonovich transformation on the partition function Z makes the fields 
$i, <£* (for charge) and £ c . appear (omitting the fields associated to S/ Ci , S c /J. We get 



with 



d$ i d$*dZ f .dt Ci Vc ia Vf ia d\ i exp - / C'(T)dT 



C'(t)=Co(t)+Ho(t)+H'j(t) 



(12) 



H'Ar) = £ ( 4, /t 



-J^£ 



Jc*iT° CT , -Js^ Ci 



fio 



A. Saddle-Point 



The saddle-point solution is obtained for space and time independent fields $o , ^o, £/o an d £ Co - In the 
magnetically-disordered regime (£/ = £ Cf) = 0), it leads to renormalized bands a and (3 as schematized in Figure 
1. Noting Oq = Jc^o*^ and = Ao, aj^lO > and fll a \0 > are the eigenstates of 

G -^(k,T)=(^+ £fe ) (13) 



<7 d T + e f 

with respectively the eigenenergies (d T + E^) and (d T + E£). In the notations: Xk = £fe — £/, J/^ = £jf — e/ 
and A fe = \fx\+ 4ctq, we get 

I^= = ( a!jk ±A fc )/2 (14) 
Let us note the matrix transforming the initial basis {c\ a , fl ) to the eigenbasis (o^ , /3iL.)- The hamiltonian 

uriuiJiM.il. Mir iu.il I l.\ (','.,, I.-. 1 1 x I. J I « i I \ . i i, T ! — i~' ' ■ — 1 '•• ••••*•••''••• — 1 

have 



being hermitian, the matrix U^a is unitary : Uk a U ka — U^U^ = 1. In the notation U\ a = I " I , we 

Ilk Vk 



G 



Uk = 



-<ro/Vk 



(15) 



i'k 



a* 



The saddle-point equations together with the conservation of the number of conduction electrons are written as 

= —Jc'YlukVk n F {E^) (16) 



O"0 



ka 



1 = W s ^ ul nF{E ^ 



Their resolution leads to 



ka 



\y F \ =De X p[-2/ (poJc)} 
Zpo°ll \Vf\ = 1 



(17) 



fi = 

where yp = M — £ f and po is the bare density of states of conduction electrons (po — 1/21? for a flat band). 
Noting y = E — ep, the density of states at the energy E is p(E) = po (l + ^o/V)- If n c < 1, the chemical 
potential is located just below the upper edge of the a-band. The system is metallic. The density of states at 
the Fermi level is strongly enhanced towards the bare density of states of conduction electrons : p(Ep)/po = 
(1 + Oo/yp) ~ l/(2po \uf\)- That corresponds to the flat part of the a-band in Figure 1. It is associated to the 
formation of a Kondo or Abrikosov-Suhl resonance pinned at the Fermi level resulting of the Kondo effect. The 
low-lying excitations are quasiparticles of large effective mass m* as observed in heavy-Fermion systems. Also 
note the presence of a hybridization gap between the a and the (3 bands. The direct gap of value 2<7o is much 
larger than the indirect gap equal to 2|j/f|- The saddle-point solution transposes to N=2 the large-N results 
obtained within the slave-boson mean- field theories [E4I-E7I. 



P 



e„=0 


2|y f || ^ 






2a 




_ r __. 


/ 




a j 









K k 



FIG. 1. Energy versus wave-vector k for the two bands a and (5. Note the presence of a direct gap of value 2<to and of an 
indirect gap of value 2 \yp\- 
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B. Gaussian fluctuations 



We now consider the gaussian fluctuations around the saddle-point solution. Following Read and Newns Q , 
we take advantage of the local U(l) gauge transformation of the lagrangian £'(r) 

$i — > Ti exp(«0i) 
fi fl exp(i9i) 
X t -»• A- + 1 00i/0r 

We use the radial gauge in which the modulus of both fields $^ and $* are the radial field r.;, and their phase ^ 
(via its time derivative) is incorporated into the Lagrange multiplier Aj which turns out to be a field. Use of the 
radial instead of the cartesian gauge bypasses the familiar complications of infrared divergences associated with 
unphysical Goldstone bosons. We let the fields fluctuate away from their saddle-point values : r, = tq + 5ri, 
Xi = Xq + SXi, £f\ = <5£ f . and £ Ci = 5£ c . . After integrating out the Grassmann variables in the partition function 
in Equation (113), we get 



Z= f V ri -DXi-D£ f .VZ c .exp[-S eff } (18) 

where the effective action is 

k,icj n i i 

with : 

\Q-lU w \] acr ' _ ( [{-i^n - fJ>)5ij ~ ty]<W - Jsififi-Taa'Sij (<7 + Jc<fri)<W Sij \ 

L "AUj ^ (o- _|_ J c 5ri)S aa <5 i: j [-iw„ + £/ + <5A,]<W<% - Jsi€ Ci - T <r<r' s ij ) 

(2) 

Expanding up to the second order in the Bose fields, one obtains the gaussian corrections S e ^j to the saddle-point 
effective action 

S Sf = \ E [( fr > ^ ) Dc 1 ^ ^) ( S ) + ( <^< ^ ) D i _1 ^ ( tfp ) 

+ ( ) D^^q,^) f g ) + ( 6£J, 6£ ) B^faw,) ( II )] (19) 

where the boson propagators split into the following charge and longitudinal spin parts 

D *(q iu)=( Jc ^~ J c(^2(q>^)+^m(q>^))] -Jc^iCq,^) \ , 2Q s 
c ^ V -JcVi(<1j*w v ) -tp ff (ti,iw v ) J 

r^W-if ■ \ I Js<p\Aq,ie}„) J s [l + Jsvlffaiw*)] 
D'i (q, iuj u ) =| . J/ n . . . „ ii . 1 



and equivalent expression for the transverse spin part D^ _1 (q, iw^). The expression of the different bubbles are 
given in the appendix. The charge boson propagator Dc(q, iuj v ) associated to the Kondo effect is equivalent 
to that obtained in the 1/N expansion theories. The originality of the approach is to simultaneously derive 
the spin propagator d| X (q, iu) v ) and Do _1 (q, iu) v ) associated to the spin fluctuation effects. Note that in the 
magnetically-disordered phase, the charge and spin contributions in S e ff are totally decoupled. 
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C. Dynamical spin susceptibility 



Next step is to consider the dynamical spin susceptibility. For that purpose, we study the linear response Mf 
to the source-term — 2S/.B (we consider B colinear to the z-axis). The effect on the partition function expressed 
in Equation Jl2] ) is to change the hamiltonian H'jij) to H'j B (t) 

iaa' \ a=x,y,z / V / i i 

(21) 

Introducing the change of variables Q — — iB/Jg, we connect the f magnetization and the ff dynamical spin 
susceptibility to the Hubbard Stratonovich fields £ * . 

1 dLnZ 
M f~~(3 dB z 

Using the expression ( p0[ ) fot the boson propagator D| _1 (q), we get for the longitudinal spin susceptibility 

X^(q, = ,, , ,, , JJ ua , (23) 



and equivalent expression for the transverse spin susceptibility xffi^^v)- The diagrammatic representation 
of Equation ( |2^ ) is reported in Figure 2. The different bubbles y//(q, iu;„), y> cc (q, iw„) and ipf c (q,iuj v ) are 
evaluated from the expressions of the Green's functions 

G//(k, iu n ) = UfeG QQ (k, iu n ) + i^G^k, iw n ) (24) 
G cc (k,iuj n ) = WfcG QQ (k, iw„) +ttfcG^(k, iw n ) 

G c /(k,iw n ) = G/ C (k,iw n ) = -u fe t; fc [G ctQ (k, kj„) - Gpp(k,i(j n )] 

where G QQ (k, iuj n ) and Gpp(\t,iuj n ) are the Green's functions associated to the eigenstates aj^lO >and /3^ CT |0 >. 
In the low frequency limit, one can easily check that the dynamical spin susceptibility may be written as 

, . s Xaa(q,«4/) + X a p(<l,iu v ) 

Xff (q,iu v ) = =2 — ; — — t — ; — r 25) 

1 - JsXaa(<h «w 1/ )x ct/3 (q, IU V ) 

for both the longitudinal and the transverse parts. 

2 n F (££) - n F (E+ +q ) 



X a/3 (q, io;„) = -= £(ufcWfc +g + wgujj 
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f c f c f c 

FIG. 2. Diagrammatic representation of Equation ( |3l| ) for the dynamical spin susceptibility x//( ( 3i 

D. Physical discussion 



From Equation (J25|), one can see that the dynamical spin susceptibility is made of two contributions 
) and Xi„ ter (q, iu u ) 

X//(q, Ui) v ) = Xintra{<\, U*>v) + Xmierfa, (26) 

with 

Xmtra(q,^) = - -5 ; — — ; r (27) 

Xi7iter(q>«*\,) = -2 t ; — ~, ; r (28) 

1 - J sX aa (q, ^)x Q/3 (q, ««m 

Xintra(q, iw„)and Xmter(q, i^>v) respectively represent the renormalized particle- hole pair excitations within the 
lower a band, and from the lower a to the upper (3 band. The latter expression is reminiscent of the behaviour 
proposed by Bernhoeft and Lonzarich [0 to explain the neutron scattering observed in UPt^ with the existence 
of both a "slow" and a "fast" component in x'(q> due to spin-orbit coupling. Also in a phenomenological 
way, the same type of feature has been suggested in the duality model developed by Kuramoto and Miyake . 
To our knowledge, the proposed approach provides the first microscopic derivation from the Kondo lattice model 
of such a behaviour. The bare intraband susceptibility x QQ (q, w) is well approximated by a lorentzian 



x a i(q>^) = Paa(q) 1 ( 1_ * r ~^y) 



(29) 



where p aa = x! (q> 0) an d To(q) is the relaxation rate of order \yp \ — Tk- This corresponds to the Lindhard 
continuum of the intraband particle-hole pair excitations x aa (<2S w) ^ as reported in Figure 3. In the same way, 
one can schematize the low- frequency behavior (u> << wo(q) of the bare interband susceptibility by 



X'^ 1 (q^)=p^(q)- 1 (l-^) CUM 



where p a p — xl^^O) and Wo(q) is a characteristic frequency-scale of the interband transitions. The value 
of wo(q) is strongly structure-dependent. In the simple case of a cubic band structure ek = —2t(cosk x + 
cos fcy+cos k z ) (tight-binding scheme including nearest-neighbor hopping), we find a weakly wavevector dependent 
frequency around q = Q of order of wo = 2|j/f| / (pqJc)- The latter result does not stand for more complicated 
band structures as obtained by de Haas- van Alphen studies combined with band structure calculations in heavy- 
Fermion compounds. In the following, we will leave o>o(q) as a parameter. Figure 3 reports the continuum of 
interband particle-hole excitations x a p" 7^ 0. Due to the presence of the hybridization gap in the density of 
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states, the latter continuum displays a gap equal to 2ao, the value of the direct gap at q = 0, and 2 \uf\, the 
value of the indirect gap at q = Q (close to kp). More precisely, we have 



W(0,w) = 4p 



UJy/u! 2 — 4(7, 



at 2<7o < us < D 



(31) 



Xafl"(Q> W ) = 2 P0 



l + ^/(2a ) : 



at 2 <u <2D 




K 2k t q 

FIG. 3. Continuum of the intra- and interband electron-hole pair excitations x aa (?> w ) 7^ and x a/j (?> w ) 7^ 0. Note the 
presence of a gap in the interband transitions equal to the indirect gap of value 2 \yp\ at q — kp, and to the direct gap of value 
2o"o at q = 0. 

Far from the antiferromagnetic wavevector Q = (tt, tt, 7t), X//(*l) w ) is dominated by the intraband transitions. 
In the low frequency limit, the frequency dependence of Xintrai^ ^ can be approximate to a lorentzian 



X//(q, w) « xLtra (q, = w 



Xmtra(q)r'mtra(q) 

+ r in i ra (q) 2 



(32) 



with 



r <ntro (q)=r (q)(l- J(q)) 



(33) 



X-intra 



Paa(q) 

(i--f(q)) 



_ /(q) = ^Xa Q (qO)x^(q^). One has: Xa Q (0,0) = p aa (0) = P {E F ) and x^O) = p . The contribu- 
tion expressed in equation ( |32| ) is consistent with the standard Fermi liquid theory. Note that the product 

1- intra 

(q)Xmtra(q) = Pace (q)r (q) is independent of I. 



Oppositely, at the antiferromagnetic wavevector Q, X//(q ; w ) is driven by the interband contribution and we 



get 



IXinter^inte 



(u - w raax ) 2 + r 2 nter 



(34) 



with 



W m ax — C0q(1 — T) 



(35) 



Winter — w o(l — -O^O 



Xznter = Pa/j/C 1 ~ I) 
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where ojq, p a p, T and I are the values of w (q), p a p{q) and T (q) and /(q) at q = Q. The role of the 
interband transitions have already been pointed out in previous works ||47| . However while the previous studies 
conclude to the presence of an inelastic peak at finite value of the frequency related to the hybridization gap 
whatever the interaction J is, we emphasize that the renormalization of x Q ^(Q,w) into Xinter{Q,,u) leads to 
a noteworthy renormalization of the interband gap. Due to the damping introduced by intraband transitions, 
Xi nter {Q,u) takes a finite value at frequency much smaller than the hybridization gap. Both the relaxation rate 
Winter vanishes and the susceptibility x'inter diverges at the antiferromagnetic transition with again the product 
Winter x'mter independent of /. Remarkably, the value w max of the maximum of Xm*er(Qi UJ )/ U} is at the same time 
pushed to zero. This excitation can be analyzed as an excitonic mode which softens at the magnetic transition. 
Such a behaviour has been effectively observed in Ce\^. x La x Ru2Si2 Jl8| with a reduction of Ti nter and w m ax 
respectively by a factor 4 and 6 when x goes from to 0.075 so when getting closer to the magnetic instability 
occuring at x = 0.08. It is likely that this mode is called to play a role in the critical phenomena observed near 
the magnetic transition. 



IV. THE SUPERSYMMETRIC APPROACH 



Traditionally, the spin is described either in fermionic or bosonic representation. If the former representation, 
used for instance in the 1/N expansion of the Anderson or the Kondo lattice models, appears to be well adapted 
for the description of the Kondo effect, it is also clear that the bosonic representation lends itself better to the 
study of local magnetism. Quite obviously the physics of heavy-Fermions is dominated by the duality between 
Kondo effect and localized moments. This constitutes the motivation to introduce a new approach to the Kondo 
lattice model (KLM) which relies on an original representation of the impurity spin 1/2 in which the different 
degrees of freedom are represented by fermionic as well as bosonic variables. The former are believed to describe 
the Fermi liquid excitations while the latter account for the residual spin degrees of freedom. 

In order to include the Fermi liquid excitations as well as the residual spin degrees of freedom, the proposition 
is to enlarge the representation of the spin operator as follows 

Sa = E to b °' + 1°' = S b + S f ( 36 ) 

aa' 

where b\ and /J are respectively bosonic and fermionic creation operators and r a (a = (+, — , z)) are Pauli ma- 
trices. Eq.(l) corresponds to a mixed fermionic-bosonic representation between Schwinger bosons and Abrikosov 
pseudo-fermions. To restrict the dimension of the Hilbert space to two, we introduce the following local con- 
straints 

n f + n b = l (37) 

The constraint restricts the Hilbert space to the two states of the form: | j) = (X&j + Yfj)\0), | J.) = 

(A"&j + F/|)|0) where X 2 + Y 2 = 1 to guarantee the state normalization to 1 and |0) represents the vacuum of 
particles: 6 CT |0) = / CT |0) = 0. X and Y are parameters controlling the weight of boson and fermion statistics 
in the representation: they will be fixed later on by the dynamics. The constraint can be viewed as a charge 
conservation of the following SU fermion-boson rotation symmetry leaving the spin operator invariant 

{f2X) = {&bl)Vl (38) 

where is an unitary supersymmetric matrix (VVt = V'V = 1). One can easily check that the representation 
satisfies the standard rules of SU(2) algebra: | f) and | J.) are eigenvectors of S 2 and S z with eigenvalues 3/4 and 
±1/2 respectively, [S ,+ ,S f_ ] = 2S Z and [S' Z ,S' ± ] = ±S' ± provided that the local constraints expressed in Eq.(2) 
are satisfied. 

In the representation introduced before, the partition function of the three-dimensional KLM can be written 
as the following path integral 
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z = E n =i 2 / Vc ia V f^VbivdvidX^ exp(- f£ dr 

With C(T) = E« ff {cld T C la + fldrfia + bldrbia) (39) 

and H = Efar e * c L c »w 

+^E i (S/ 4 + S 6j ).s i -ME^c 4 

Note the presence of two terms in Z coming from the contributions of the states satisfying Q$ = Sf + (— l) n /2 = 
respectively for n = 1 and 2. The time- independent Lagrange multipliers X\ and ^ are introduced to enforce 

(n) 

the local constraints Q\ = and Pj = + — 1 = 0. Performing a Hubbard-Stratonovich transformation 
and neglecting the space and time dependence of the fields in a self-consistent saddle-point approximation, we 
have 

z = En=i,2 fdr]dr]*C n (<r , A , ?7, V*) z n(v, V*) 
Zniv, V*) = E CT / Vc ia Vf ia Vb ia exp (- jj dr(£(r) + W'„ CT ) 



with w / w = e*(/L4.*L)^ 



j/ + aA n) /2 <r 





(40) 



where C n (ao, Ao, 77, 77*) is an integration constant. £/ and Ag"^ are the saddle-point values of the Lagrange 

(n) 

multipliers vi and \\ . Note the presence of a Grassmannian coupling r\ between a a and in addition to the 
usual coupling ao between Cj CT and responsible for the Kondo effect. In the following, Hq is indifferently used 

for any Hq°. Hq is of the type ( ^ ^ ^ in which a, b (p,a) are matrices consisting of commuting (anticommuting) 

variables. Note the supersymmetric structure of the matrix Hq similar to the supermatrices appearing in the 
theory of disordered metals . 

Hq being hermitian, the matrix U< transforming the original basis ijft — (/t, c' , b') to the basis of eigenvectors 
& = (a^,/3^,7^) is unitary (UW — U^U — 1). & — ip^U< with U< a supersymmetric matrix, ofi and (ft are 
the fermionic eigenvectors whose eigenvalues, determined from det [(a — E) — a(b — E)^ 1 p\ = 0, are 



T 2 

7^ is the bosonic eigenvector whose eigenvalue, determined from det [(b — E) — p(a — E)~ 1 a~\ = is E 7 — ef. 

In this scheme, ao and Ao are slow variables that we determine by solving saddle-point equations, while 77, 77* 
are fast variables defined by a local approximation. As we will see, the latter approximation incorporates part 
of the fluctuation effects. Indeed, performing the functional integration of Eq.(J40|) over the fermion and boson 
fields |Q yields a superdeterminant (SDet) form written as follows 

Z( V ,r 1 *) = SDet(d T + H ) , 

™ ^ Q rr, DetiG^-aDp) (41) 
where SDet(0 T + H)= 1__ ^ , 

G" 1 = <9 r + a and D" 1 = <9 T + 6 . 

Expanding to second order in 77, 77* allows us to define the propagator G^* (k, iu) n ) associated to the Grassmann 
variable 77 and hence the closure relation for Xq = (7777*) 



G m *{k,icj n ) , (42) 
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with ^(k,^)= [1 _ jn / (k)ia;n) 



and U° cb = - ^2 G cc(k + q, iiv n + icj u )D(q, iuj n ) . 



Contrary to |49| which assumes x\ = leading to a two-fluid model description, the closure equation Eq.(42) 
defines a finite x^. This parameter x^ plays a major role in controlling the relative weights of fcrmion and boson 
statistics. It is directly connected to the X and Y parameters introduced in the initial representation of the states 
: A 2 = xl/{al + x 2 ) and F 2 = og/(og + *§). 

The resolution of the saddle-point equations, keeping the number of particles conserved, leads to 

y F = -Uexp[-l/(2Jpo)] , 
1 = 2p (g 2 

2n ' (43) 
A« = A^= ,' 

where j/f = fi — £f and po = 1/22? is the bare density of states of conduction electrons. From that set of 
equations, we find : £/ = 0. 

The resulting spectrum of energies is schematized in Figure |ij At zero temperature, only the lowest band 
a is filled with an enhancement of the density of states at the Fermi level (and hence of the mass) unchanged 

from the standard slave-boson theories P^ F ^ = 1 + (^o + x o) _ ^ _j_ i This large mass enhancement 

Po Vf _ K-Vf) 

is related to the flat part of the a band associated with the formation of the Abrikosov-Suhl resonance pinned 

at the Fermi level. While this feature was already present in the purely fermionic description, it is to be noted 
that the formation of a dispcrsionless bosonic band within the hybridization gap is an entirely new result of the 
theory. 




k F k 

FIG. 4. Sketch of energy versus wave number k for the three bands a, /3, 7 resulting of the diagonalization of supersymmetric 
H . 



The relative weight of boson and fermion statistics in the spin representation is related to Xq : rib/nf — Xq/(Tq. 
It is then interesting to follow the J-dependence of Xq as determined by the closure equation (|42]). The result is 
reported in Fig. ^. This bell-shaped curve can be interpreted in the light of the exhaustion principle mentioned 
in the introduction. In the limit of large J, the Kondo temperature-scale Tjc = Dexp[— 1/(2 Jpo)] is of order of 
the bandwidth. One then expects a complete Kondo screening as can be checked by remarking that the weight 
of c in the a quasiparticle at the Fermi level (noted vf. F ) just equals the added weights of f and b at the Fermi 
level (respectively noted u| and p\): v1 f /{u\ f + p\) = y F /{<Jo + x%) = 1. The Kondo effect being complete in 
that limit, there is no residual unscreened spin degrees of freedom: it is then natural to derive a zero value of Xq 
(and hence of rib). The opposite limit at small J corresponds to the free case of uncoupled impurity spins and 
conduction electrons. It also leads to: Xq — 0. The finite value of Xq between these two limits with a maximum 
reflects the incomplete Kondo screening effect in the Kondo lattice, the unscreened spin degrees of freedom being 
described by bosons. 
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FIG. 5. J/-D-dependence of the coupling x 2 , — {r)rf) fixing the relative weight of fermion and boson statistics. The unit on 
the vertical-axis is D 2 . 

Largely discussed in the litteraturc po|] is the question concerning the Fermi surface sum rule: do the localized 
spins of the Kondo lattice contribute to the counting of states within the Fermi surface or do they not? Depending 
on the answer, one expects large or small Fermi surfaces. The supersymmetric theory leads to a firm conclusion 
in favour of the former. One can check that the number of states within the Fermi surface is just equal to 
n c + Tit, + rif, i.e. n c + 1. The Fermi surface volume includes a contribution of one atate per localized spin in 
addition to that of conduction electrons (2(].|5(|. The latter conclusion appears sensible if one recalls that the 
KLM is an effective hamiltonian derived from the periodic Anderson model (PAM). 

Let us now consider the response functions to some external fields namely the dynamical spin suscep- 
tibility x afc (q, oj) and the frequency-dependent optical conductivity a ab (ui) (a,b — x,y,z). For that pur- 
pose, we introduce the Matsubara correlation functions associated with the operator O a (q, r): x ab (Qj*^) = 
J^dr exp^" T (T T O Q (q,r)O b (-q,0)). The operator related to the spin-spin correlation function is the 
a-component of the spin expressed in the mixed representation introduced in the paper by S* a (q) = 
Sfc a a 1 fl+q a T a<j' fk.<r' + &fc+ g cr T acr' bk,a' • As usua l) the dynamical spin susceptibility is then derived from the 
spin-spin correlation function by the analytical continuation iu> u — ► lu + i0 + . In the same way, the opera- 
tor related to the current-current correlation function is the a-component of the c-current. In the case of 
a cubic lattice: J*(q) = 2 a a , sin k a c ! k+q a Ck.a- The frequency-dependent optical conductivity is then 
obtained from the current-current correlation function by the analytical continuation following: a ab (ui) = 
[x ab (q,w + i0+)-x a& (q,*0+)] /ico. 





0.00 0,02 0,04 0,06 0,08 0,10 

IB 

FIG. 6. Frequency-dependence of the optical conductivity u(lj) at T = for D = 0.8 and Tk = 0.001. 

By expanding the previous expressions in the basis of the eigenstates (a^tyn f we have computed the 
frequency dependence of x a& (Q, at the antiferromagnetic wavevector Q and a ab (uj) at zero temperature. The 
two response functions show very different frequency dependence. The frequency-scale at which the dynamical 
spin susceptibility takes noticeable values is much smaller than for the optical conductivity. This can be under- 
stood in the following way. The bosonic 7-band is called to play a role only when spin is concerned namely for 
the dynamical spin susceptibility. That feature comes from the fact that the spin is related to both fermionic 
and bosonic operators while the c-current is simply expressed within fermionic operators. Therefore, one can 
show that the dynamical spin susceptibility involves transitions between all three bands a, (3 and 7. The main 
contribution for x ab (Q J ti; ) is due to the particle-hole pair excitations from the fermionic a to the bosonic 7 
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band. Oppositely, the optical conductivity is associated with transitions between fermionic bands only. As can 
be seen in Figure ||, a gap appears in the frequency dependence of cr(w) equal to the direct gap between the a 
and (3 bands. The latter result agrees with the predictions of the dynamical mean-field theory in the limit of 
infinite dimensions f3l| ] . The whole discussion clarifies the physical content of the novel bosonic mode brought 
by the supersymmetric approach. That mode is related to the spin excitations. It introduces new features in 
the dynamical spin susceptibility by comparison to the standard slave-boson theories while it does not affect the 
optical conductivity. 



V. CONCLUSION 



Important progress has been made the last years in the understanding of the Kondo lattice model with the 
development of new functional integral approaches. They have enlightened as to the nature of the ground state 
and the existence of collective modes. They open up new prospects for the description of the critical phenomena 
associated to the quantum phase transition in Heavy-Fermion systems. A complete study of the quantum phase 
transition will probably requires the use of the Group Renormalization techniques for which the functional 
integral approaches presented here might constitute the framework. 
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APPENDIX 

The expressions of the different bubbles appearing in the expression of the boson propagators (cf. Eq.po|) are 
given here (with i=l, 2, m or ff) 

^i(q, iu v ) = ipi(q, iuj v ) + <Pi(—q., —iu) v ) (44) 
<pi(q, iu v ) 

</? 2 (q, ii» v ) 

¥> m (q, iu) v ) 



tp l f c (q, iuJ„) 



= -pYl G c/ ( k + q.^n + iw v )G'} fo Qs,iu n ) 
= ~\ H G£jk + q>„ + ^)G? /o (k,w n ) 

k(T,iuj n 
ka,iuj n 

= -4 G'ff Os- + q.,iuJ n + itJ u )Gf fo (k,iu) n ) 

kcr : iuj n 
kcr,iuj n 

= —S G /c ( k + q.^n + iuj„)G a fCQ (k,iuj n ) 



1G 



ipj f (q,iuj u ) = -- ^ G/ /o (k + q,^„ + iw l/ )G} /o (k,iw„) 

k(j,iuj n 

(p^.(q,iuj v ) = -- ^ Gl C0 (k + q,iu n + iu„)G l CCo (k,iuj n ) 

kcr,iuj n 

(pf c (q,iuj v ) = -- ^ G^ Co {k + q,iu} n + iuj u )G l fco (k,iuj n ) 

where G^ Co (k, icu n ), G^ g (k,ioJn) and Gy C[) (k, iuj n ) are the Green's functions at the saddle-point level obtained 
by inversing the matrix Gg(k, r) defined in Equation (|l3|). 
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